Energy resolved STM mapping of Ceo on metal surfaces: A theoretical study 
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We present a detailed theoretical study of scanning tunneling imaging and spectroscopy of Ceo on 
silver and gold surfaces, motivated by the recent experiments and discussion by X. Lu et al. [PRL 
90,096802 (2003) and PRB 70, 115418 (2004)]. The surface/sample/tip system is described within 
a self-consistent DFT based tight-binding model. The topographic and conductance images are 
computed at constant current from a full self-consistent transport theory based on nonequilibrium 
Green's functions and compared with those simulated from the local density of states. The molecular 
orbitals of Ceo are clearly identified in the energy resolved maps, in close correspondence with the 
experimental results. We show how the tip structure and orientation can affect the images. In 
particular, we consider the effects of truncated tips on the energy resolved maps. 



PACS numbers: 73.63.-b,73.22.-f,73.40.Gk 

I. INTRODUCTION 

Scanning tunneling microscopy (STM) was introduced 
to image surface structures in real spacei, and is nowa- 
days used to obtain more complete information on the 
local electronic properties with the scanning tunneling 
spectroscopy (STS) technique. Its operation is based on 
the measurement of the tiny tunneling current flowing 
from a very sharp conducting tip to the surface when a 
bias voltage is applied. Images are obtained by record- 
ing the vertical displacement of the tip as it is moved 
across the surface, while keeping the tunneling current 
constant with a precise feedback mechanism. If the cur- 
rent increases due to a protrusion on the surface, the tip 
is lifted to bring the current back to the set value, there- 
fore revealing the topography of the substrate. 

Tersoff and Hamann 2 have shown how the images do 
not only reflect the geometric structure of the surface, 
but also depend on the electronic density of states of 
the sample, which is identified, to a first approximation, 
with the differential conductance. This observation has 
opened the way to spectroscopic measurements with a 
spatial resolution which is by far not accessible by other 
surface science techniques. This capability is of particu- 
lar interest in the study of adsorbates at surfaces, which 
can be addressed individually by the probe. The STS 
signal, on the other hand, is strongly influenced by the 
nature of the contact and the structure of the tip, and the 
extraction of information on the physical and chemical 
properties of the measured sample is possible only with 
the support of electronic structure calculations with an 
atomistic description of the system 3 ^^. 

In this paper we discuss from a theoretical point of 
view the powerful combination between microscopy and 
spectroscopy which is peculiar of STM/STS. The ac- 



curacy of recent STM experiments on Cgo monomers 
on metallic surfaces^ represents a rigorous benchmark 
for any theoretical modeling of STM imaging and spec- 
troscopy. Quite interestingly, energy resolved mapping 
of the molecular orbitals of adsorbed Ceo monomers has 
already been used as a control reference for even more so- 
phisticated experiments 8 ; similarly, we would like to dis- 
cuss in detail STM/STS mapping of this system based on 
Green's function methods, in view of the investigation of 
adsorbed molecules and clusters- 9 which have not yet been 
experimentally characterized with this technique. Topo- 
graphic and energy resolved maps are computed with a 
realistic model of the STM probe, with the simulation 
of the tip trajectory over constant current surfaces; the 
effects of the tip structure on the rendered images are 
illustrated in detail, showing that truncated tips may du- 
plicate the characteristic features depending on the ori- 
entation of the topmost atoms with respect to the sam- 
ple. Also, we refine previous calculations based solely 
on the local density of states (LDOS) of the supported 
molecule^. We confirm that the images obtained from 
the LDOS can effectively approximate the patterns ob- 
served in the STS maps, with the additional piece of ad- 
vice that the simulated probe trajectory should be taken 
at the correct distance from the molecule. 

The paper is organized as follows. The theoretical 
framework is briefly reviewed in Sec.[n] The experienced 
reader may skip this part and move directly to the re- 
sults in Sec . IIIII Equilibrium electronic properties of the 
supported Cgo are discussed in Sec. lIII Al while the cal- 
culations based on transport theory are in Sec. lIII Bl ex- 
tended by simulation of topographic and spectroscopic 
images with different types of tips in Sec. lIII 01 Finally, 
in Sec. lIII Dl we specify how STS maps can be modeled 
using the local density of states. 
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II. THEORY 

We consider a fully atomistic description of surface, 
molecule and STM probe. The geometry of substrate 
and the STM electrode are fixed from bulk properties, i.e. 
surface relaxation effects are neglected, and the structure 
of the molecule is assumed to remain unchanged after 
adsorption. The transport properties are computed with 
a microscopic approach using non-equilibrium Green's 
functions, and the electronic structure is modeled within 
a single-particle picture, based on a self-consistent tight- 
binding (TB) model which is parameterized from density 
functional theory (DFT). 



A. TB— DFT model 

Electron wavefunctions are approximated by an ex- 
pansion in an atomic-like orbital basis^, which is de- 
termined by solving the Kohn-Sham equations of the in- 
dividual atoms within the local density approximation 
(LDA). The TB Hamilton matrix expressed in the local- 
ized basis consists only of two-center integrals which can 
be obtained exactly by using a Slater type functions for 
the atomic orbitals. This scheme is refined to include 
charge self-consistency by expanding the DFT energy 
functional to the second order in the charge-density fluc- 
tuationsii*i£. Such an expansion leads to the appearance 
of an additional term to complete the TB Hamiltonian 
that contains the contribution of the net atomic charges 
on the local effective potential. Like in a full DFT calcu- 
lation, self-consistency is obtained with iterative meth- 
ods, such that the local density, given by the sum of 
the (approximated) Kohn-Sham eigenfunctions, is con- 
sistent with the effective potential. The main advantage 
of this approach lays in its simplicity and computational 
efficiency, while basic key effects like charge transfer and 
redistribution are still handled correctly. 



B. Equilibrium electronic properties of adsorbates 

Inglesfieldiiiii introduced the concept of embedding 
potential to perform real space electronic structure cal- 
culations in a domain surrounding impurities or defects in 
an otherwise perfect crystal. The form of the embedding 
potential is derived from the Green's function of the solid, 
and is included in an effective Schrodinger equation, to be 
solved only in the region of interest, ensuring the correct 
matching of the solutions with the wavefunction in the 
rest of the solid. A complementary approach considers 
the defect as a perturbation of the ideal system, and has 
been first proposed by Williams and Langi^. A localized 
atomic orbital basis {(^(j))} * s introduced to represent 
in matrix form the Hamiltonian H and the (retarded) 
Green's function G, defined by 



S = Sfj, v — {tp^\tp v ) being the overlap matrix. The indices 
for matrices with a tilde run over the atomic orbitals of 
the surface and the cluster. 

Measurable quantities can be computed directly from 
the Green's function; for example, from the spectral den- 
sity matrix 



;Im [Gp V (E)\ 



(2) 



(the factor 2 is explicity included to point out that we 
consider a spin-unpolarized system), it is possible to 
compute the local density of states 



and therefore the density 



/E-p 
dE p(r,E); 
-OO 



(3) 



(4) 



where E-p is the Fermi level (degenerate case). 

The localized basis allows to translate the cluster- 
surface partitioning of the system into the block structure 
for the overlap, the Hamilton, and the Green's function 
matrices of Eq. . The Green's function of the isolated 
cluster G° is modified by the presence of the surface, and 
can be obtained from the solution of the Dyson equation 



G = G° + G°SG 



(5) 



where G refers to the Green's function of the "embed- 
ded" cluster, and £ is the matrix analogue^ of Ingles- 
field's embedding potential. Note that, in contrast to 
Eq. 0J, all matrices in Eq. are defined for orbitals at 
the cluster only. X can be expressed as a self-energy 



Z(E) = V(E)G° S (E)VHE), 



(6) 



(E+irj) S — H G(E) = 1, r\ -> H 



(1) 



thereby including the scattering processes due to the in- 
teraction with the surface, given by the combination of 
the unperturbed surface Green's function G° and the 
tunneling matrix V = E S cs — H cs . S cs and H cs are 
the overlap and Hamilton matrices between orbitals at 
the surface and the cluster. For an orthogonal basis S cs 
would vanish and V would be independent of the energy 
E. 

Having chosen a localized orbital basis, only a limited 
number of elements in the tunneling matrix are different 
from zero, and Gj? needs to be evaluated just for the sur- 
face atoms which are within a finite distance from the 
supported cluster. The two-dimensional periodic struc- 
ture of the surface allows to reduce the corresponding 
overlap and Hamilton matrices in a block-diagonal form, 
based on the transformation to the two-dimensional mo- 
mentum space. The surface is then partitioned in a series 
of adjacent layers, and the Green's function for the first 
layer, the needed submatrix of G° , can be efficiently com- 
puted for any k-point using recursive methods^. The 
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conversion to the real space representation is finally ob- 
tained by integrating over a finite number of k-values in 
the first Brillouin zone (we sampled over a regular grid 
with at least 1600 points). 

Charge self-consistency requires that the effective po- 
tential entering in the system hamiltonian leads to a 
charge density which generates the same effective poten- 
tial. It is obtained with an iterative scheme, by com- 
puting the local charge densities from G, to be fed into 
the local effective potential of the cluster Hamiltonian 12 , 
therefore leading to a new cluster Green's function, which 
provides the density for the next iteration, until conver- 
gence is reached. The integral of Eq. is evaluated 
most efficiently taking advantage of the analytic struc- 
ture of the Green's function, with a contour integration 
in the upper complex pian o 15 ! 18 . Since the substrate 
is described as an ideal surface, the effect of screening 
due to the charge redistribution in the metal is achieved 
by including image charges, which generate an external 
field acting on the cluster region, and provide the correct 
boundary conditions for the Poisson equation, such that 
the electrostatic potential on the first atomic plane of the 
surface is uniform and set to zero. 



C. Transport 

The expression for the density in Eq. is valid only 
in equilibrium conditions. When the sample is probed 
with an STM tip at nonzero bias V, the local charge is 
obtained from the correlation function G < defined byi^ 

G< = GE^, (7) 

where the lesser self-energy 

s< = i[/(^- Ms )r s + /(^ Mt )r t ] , (8) 

depends on the chemical potentials [x s = Ep and /it = 
fi a + \e\V = Ep + \e\V of the surface and the tip, respec- 
tively. In Eq. JSJ, / is the Fermi function, and T s and T* 
are the broadening matrices. They include the interac- 
tion with the surface and the tip through the respective 
self-energies 

r s = i [s s - [s 8 ] 1 '] , r* = i [e* - [s 4 ] 1 ] , (9) 

where the cluster-tip self-energy, £*, is defined as for the 
surface case in Eq. © ■ The non-equilibrium analogue of 
Eq. J2J reads 

i r°° 

n(r, V) = - dE G< U (E, V)^(t)^(t), (10) 

and enters in self-consistent relaxation of the electronic 
charge in the cluster, as for the equilibrium case. The 
correlation function, and therefore the non-equilibrium 
density, depends directly on the applied bias V due to 
the shift of the tip chemical potential /i t . The numerical 
procedure for the computation of the electron density 



splits the integral in Eq. I|10|) into equilibrium and non- 
equilibrium contributions, exploiting the same contour 
integration used to evaluate the equilibrium density 18 . 

In the crudest approximation, the tip electrode is as- 
sumed to be a flat metal surface, to be treated exactly 
like the one supporting the sample. The boundary con- 
ditions for the Poisson equation are that the surfaces of 
the two electrodes have a constant uniform potential dif- 
fering by the bias voltage, and the correct electrostatics 
is obtained using a series of image charges on both sides. 
This choice of boundary conditions implies that G< de- 
pends on the external potential generated by the two 
electrodes, which enters both in the cluster and tunnel- 
ing Hamiltonians, and therefore determines the cluster 
Green's function, now G(E, V), and the broadening ma- 
trices. Realistic tips are easily included into the model 
by considering the tip atoms as part of the cluster, which 
will then just be a larger sample sandwiched between the 
two parallel electrodes. The arrangement of the tip atoms 
was assumed to be consistent with the crystal structure 
of the underlying ideal surface. 

In the case of coherent transport, i.e. when there is 
no inelastic scattering within the sample, the stationary 
current is given by the simple formula^ 

2e f°° 

I = J dE T(E, V) [fiE-fH) - f(E-^ s )} (11) 

where T is the transmission function 21 

T(E, V) = Tr [G(E, V)T S (E)G^(E, V^iE-^V)] . 

(12) 

The effects of scattering processes within the conductor 
due to electron-electron or electron-phonon interactions 
can be treated consistently within the non-equilibrium 
Green's function framework by adding the correspond- 
ing self-energy termsi2iS£, but will not be considered in 
this paper. Electron-electron interactions are incorpo- 
rated here at a mean-field level, which does not cause 
incoherent effects^. 



III. RESULTS 

For the calculations presented below, the Cgo was 
placed with different orientations on Ag(lll) and 
Au(lll) surfaces at a distance of 2.4A, consistent with 
the value obtained from first principles calculations for 
monolayers^. The structure of the supported molecule 
was not relaxed from the free configuration. 



A. Total DOS of supported cluster 

Figure^] shows the total DOS of Cqq for three different 
orientations on Au(lll) and Ag(lll) surfaces. For each 
of these orientations it was observed that the number and 
position of the peaks remain unchanged as the molecule 
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FIG. 1: Total DOS p(E) of Ceo with different orientations on Au(lll), left column, and Ag(lll), right column. The top figure 
shows the position of the HOMO, LUMO and LUMO+f orbitals for the free molecule relative to the Fermi level of the Au and 
the Ag surface, respectively. 




FIG. 2: Representation of the junction formed by the surface/Ceo /tip/electrode, a: flat electrode (no tip atom), b: sharp tip 
(single atom at tip end), c: truncated tip (three atoms at tip end). 



is translated in the plane parallel to the surface, show- 
ing that in both cases there is no chemical bonding with 
the substrate. For C 60 on Ag(lll), besides the broad- 
ening and splitting of the levels of the free molecule due 
to the interaction with the surface, we observe a shift of 
+0.7 eV associated with a net transfer of 0.3 electrons 
to the molecule, as given by the Mulliken charge analy- 
sis. In the case of the Au substrate there is almost no 
charge transfer (0.01 electrons), which leaves the posi- 
tion of the free cluster energy levels unchanged. Also, 
we observe that the LUMO peak for the Ag case appears 
much more broadened than for Au, with its lower en- 
ergy tail approaching the Fermi level, similarly to what 
was observed by other DFT calculations&t The energy 
difference -Elumo ~ -Ehomo is equal to 2.0 eV in both 



cases, while the experimental gap is 2.0 eV for Ag(100) 
and 2.7 eV for Au(lll)££ The discrepancy is to be im- 
puted mainly to the intrinsic limits of LDA eigenvalues 
as a representation of electronic excitation energies. 



B. Transport calculations 

For the calculation of the conductance as a function 
of energy and of the current -voltage characteristics we 
consider the combination Cgn/tip as a single conduc- 
tor placed between two parallel, infinite, ideal surfaces 
(Fig. |2| . This way of partitioning the system is conve- 
nient since the Green's function of the surfaces can be 
computed efficiently (see Sec. lII Bll . We choose Au as el- 
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FIG. 3: Comparison of spectral density, conductance and differential conductance for Ceo on Au(lll) and Ag(lll) surfaces. 
Bottom panels: Total DOS p(E). Center and top panels: Conductance G (dashed lines) and differential conductance dl/dV 
(solid lines) for a probe made with the flat electrode (Fig.Hk) or the sharp tip (Fig.[2J>) , respectively, each placed at 5.5 A from 
the top of the molecule. Note the different scales for the peaks of G and dl/dV for E — Ef < —0.5 eV and E — Ef > —0.5 eV, 
respectively (Go = 2e 2 /h). 




FIG. 4: Computed topographic images of Ceo on Au(lll) with different orientations; the topmost features are: (a) a hexagon 
ring, (b) a 6:6 bond, (c) a 5:6 bond, (d) an apex atom, and (e) a pentagon ring. Current was set to 77.5 pA, the bias to +1.8 V. 
The size of each frame is 20 x 20 A. Compare to Fig. 1 in Lu et alA 



ement for our probe, which, in contrast to the transition 
metals used in experiments, has a fairly uniform DOS 
near the Fermi level, without any contribution from the 
d states. This choice allows us to single out the effects 



of the tip atomic structure on the calculated maps. We 
will use the term conductance for the quantity 

G(E) = G T(E,V=0), (13) 
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with G = 2e 2 /h and T(E,V) defined in Eq. (O; the 
transmission function is therefore computed from the 
self-consistent Hamiltonian of the system at zero bias. In 
contrast, the differential conductance dl/dV is obtained 
from the derivative of the total current (fill) for each bias 
value, and will be sensitive to the external field, leading 
most notably to shifts of the energy levels and therefore 
to the peak positions. 

In Fig.El we compare the spectral density p{E) of the 
supported molecule with the conductance G and differ- 
ential conductance dl /dV curves as a function of the 
bias. The latter two are computed both for a flat elec- 
trode (Fig.[2k) and a sharp tip consisting of ten atoms 
(Fig. I33). Conductance and total DOS peaks appear at 
the same position. However, the relative height of the 
peaks changes dramatically; when switching to the sharp 
tip, the conduction through the HOMO level drops by 
about an order of magnitude, while for the LUMO and 
LUMO+1 there is a slight increase in the conductance. 

At finite bias, the external field modifies the local 
charge distribution and displaces the peaks of the differ- 
ential conductance with respect to the equilibrium energy 
levels. Due to the large separation between the tip and 
the Ceo, these shifts can be appreciated only at large bi- 
ases. A similar behavior is observed for a sharp tip, which 
in this case was placed centrally above the Ceo molecule. 
Here, as well, the HOMO conducts current better than 
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FIG. 5: Dependence of calculated STS spectra on the tip 
position for the Ceo at the Ag(lll) surface with a 6:6 bond 
on the top. The spectra are shifted along the vertical axis for 
clarity. 



the LUMO and LUMO+1 states; the difference in con- 
ductance, however, is attenuated. This behavior depends 
in fact on the position of the tip, and will be discussed 
in the following section. 



C. STM/STS imaging with realistic tip models 

1, Normal maps 

We first model STM/STS imaging with a sharp tip 
(Fig. \2p). The tip is moved on a three-dimensional 
finite-element mesh placed above the molecule, with a 
grid spacing of 0.5 A. For each position the current I 
and the differential conductance dl/dV are computed 
using Eq. (|11(1 . The constant -current iso-surface is ob- 
tained by interpolation, and provides the topographic 
image of the molecule at the given voltage bias, as shown 
in Fig.QJi-e for five different orientations. The bias is 
set to +1.8 V, and the pentagon rings of the C 60 bulge 
out in a doughnut-like shape from the spherical structure 
of the molecule, allowing the direct identification of its 
orientation on the surface. 




FIG. 6: Energy resolved maps of molecular orbitals for Ceo 
on Au(lll) and Ag(lll) surfaces with 6:6 bond as the up- 
permost feature. The maps for HOMO, LUMOa, LUMO/?, 
and LUMO+1 are obtained at -1.32 eV, 0.40 eV, 0.64 eV and 
1.46 eV for Au(lll) and at -1.40 eV, 0.30 eV, 0.55 eV and 
1.36 eV for Ag(lll). The size of each frame is 20 x 20 A. 
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Based on the topographic scan of Ceo with a 6:6 bond 
(a carbon bond shared by two hexagon rings) as a top- 
most feature reported in Fig.^fc, we show in Fig.El con- 
duction curves for six tip positions. Heights and shapes 
of the peaks vary considerably as the tip is moved over 
the molecule, and we also observe slight shifts in their 
position in the energy scale, induced by the changes in 
the local electric field surrounding the STM probe. 

The strong spatial dependence of the conductance al- 
lows to construct energy-resolved maps (Fig.EJ to rep- 
resent the states associated with each peak. The sim- 
ulated STS maps are rendered by displaying the com- 
puted dl/dV values with a contour plot on the three- 
dimensional tip trajectory, which is then projected on the 
two dimensional plane of view, which in this paper will 
always be parallel to the metal surface. The grayscale 
code associates bright areas to high conductance, dark 
areas to low conductance. 

The structure of the HOMO, LUMO and LUMO+1 
molecular orbitals is clearly visible in the calculated maps 
reported in Fig-EI similarly to what was observed in ex- 
periments (cf. Fig. 3 of Lu et al.— and Fig. 3 of Lu et al.— ). 
In particular, the LUMO level appears to be composed 



of two orbitals, LUMOcv and LUMO/3, having a comple- 
mentary spatial distribution of the conductance. This 
pattern was observed in experiments 6 for the Ag(100) 
surface, where the LUMO peak displays a marked split- 
ting of 0.4 eV. By calculating maps for two energies (0.40 
and 0.64 eV) within the LUMO peak, we are able to ex- 
tract two conjugate LUMO images also for the case of the 
Au substrate, which points out that the two states sim- 
ply arise from the lifting of the degeneracy of the LUMO 
level of the free molecule due to the contact with the sur- 
face. We therefore conclude that the splitting was not 
observable in the case of the Au surface^ due to the lim- 
its in resolution for STM spectroscopy. What remains 
still unclear is the reason why the states are so well sep- 
arated in the experiments performed with the Ag(100) 
substrate 6 . The main clue is the large charge transfer 
to the molecule, compared to the neutrality for the case 
of the Au surface, which may induce a deformation in 
the molecule therefore lifting degeneracy of the LUMO 
level 6 .. 



2. Abormal maps 
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FIG. 7: Energy resolved maps with a truncated tip of Ceo on 
Au(lll) with the 6:6 bond on top (compare Fig.HjJl. Column 
a): tip in the reference position; column b): tip rotated by 
23° degrees around the axis normal to the surface. The ori- 
entation of the three atoms at the bottom of the STM probe 
is reported for each series of images. The size of each frame 
is 20 x 20 A. 



In order to investigate the sensitivity of the energy re- 
solved maps to the structure of the probe, we considered 
a truncated tip, having three atoms at its end (Fig.|2t). 
Figure[7| shows the calculated images for Ceo with two 
different orientations of the tip (columns a and b), whose 
position is given by the center of mass of the three end 
atoms. At all energies, the overall shape of the maps 
reflects the triangular arrangement of the three atoms 
at the end of the truncated tip; a rotation of the probe 
around the axis normal to the surface rotates the image 
of the molecule (compare columns a and b of Fig.0) , de- 
spite the fact that the position of the molecule remains 
unchanged. For the LUMO+1 level we observe that the 
bright rings corresponding to the pentagons in the Cqq 
are duplicated, which should be linked with the presence 
of two adjacent atoms in the triangle. For the orienta- 
tion a) the alignment of the tip is such that this sort of 
dichroism is present only for the left side of the maps, 
while when the tip moves to the right side only one atom 
remains closer to the molecule, and the single bright ring 
is recovered. A similar behavior can be noticed also for 
LUMOa and the complementary map LUMO/3. 

Most features described so far can be roughly explained 
assuming that, like in the case of the sharp tip, at every 
position of the simulated scan tunneling occurs between 
the sample and only one atom of the probe, namely the 
one closest to the adsorbed molecule. The main differ- 
ence is that for the truncated tip such atom is not always 
the same, which leads to the dichroism mentioned above. 
Only in the central region all three atoms may contribute 
equally; the bright spot appearing right at the center of 
the LUMO+1 map in Fig.0 appears precisely because of 
the rise in conductance given by the sum of the tunnel- 
ing currents through the three end atoms. We also note 
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FIG. 8: Energy resolved maps from the LDOS for Ceo on 
Ag(lll) with a 6:6 bond on top. The three columns from 
left to right show simulated tip trajectories which are approx- 
imately 2, 3.5, and 5 A away from the Ceo- The size of each 
frame is 20 x 20 A. 



that, compared to the well contrasted images of Fig. [51 
the coherent structure of the HOMO level becomes less 
evident when imaging with the blunt tip. In this case, 
the coarser nature of the probe sensibly changes the tra- 
jectory of the tip, which therefore moves over a surface 
having a lower corrugation, borrowing form the terminol- 
ogy used for fiat surfaces. The tunneling current involves 
more than one atom of the tip, hindering the resolution 
of the fine structure of the HOMO orbital. 

These results clearly demonstrate how the details of 
the energy resolved maps of Ceo can reveal the structure 
and orientation of the STM tip, and how this kind of 
analysis could serve as a tool to characterize the proper- 
ties of the probes used in experiments. 



Equation (|14fl allows to calculate the tip trajectory, based 
on the assumption that the convolution effects due to the 
tip shape and states can be neglected. The local differen- 
tial conductance d/(r, y)/dV tx p(r, E) can be projected 
onto the iso-current surface to simulate the experimental 
procedure followed to measure STS maps. The frames 
in Fig.|Hl show that the patterns in the molecular or- 
bitals discussed so far can be reproduced with the LDOS, 
provided that the calculated tip trajectory is placed far 
enough from the molecule. The standard choice is in fact 
to consider the LDOS at a few A from the atoms of the 
sample, and we verified that the peculiar LUMOa and 
LUMO/? images can be obtained only considering iso- 
surfaces at 3.5-5 A from the Cgo- 



IV. SUMMARY 

We have presented a detailed numerical investigation 
of STM topographic and energy-resolved mappings of 
Cgo on gold and silver surfaces. The transport calcu- 
lations are based on non-equilibrium Green's function 
methods, with a DFT-based tight-binding model. Self- 
consistency allows to study charge redistribution on the 
molecule due to the presence of the supporting surface, 
as well as due to external field generated by the voltage 
difference between the substrate and the STM probe. 

We have discussed the importance of realistic models 
for the tip; using a sharp tip we could perfectly reproduce 
recently measured energy resolved maps for Ceo on metal 
surfaces 6 - 7 . The images are obtained by projecting the 
contour plot of the differential conductance mapped on 
the three-dimensional surface described by the tip as it 
scans the sample with a fixed tunneling current. Images 
of a similar quality can be obtained from the LDOS using 
Tersoff-Hamann theory to calculate the conductance if 
the simulated constant-current surface is sufficiently far 
from the cluster (here more than 3 A). 

Additionally, we have modeled the STM/STS mapping 
with a truncated tip, and described the appearance of ab- 
normal images that allow to reconstruct the structure of 
the probe by direct comparison with the ideal maps. We 
finally suggest that this procedure could be used to gauge 
the properties of the STM probes, which is generally un- 
known in STM experiments. 



D. STM/STS imaging from LDOS 

Within Tersoff-Hamann theory^, the STM topo- 
graphic and spectral images are modeled assuming that 
the current at bias V is proportional to the spatially re- 
solved density of states: 

pE F + \e\V 

I{r,V) oc / dE p{r,E). (14) 

Je f 
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